Ferromagnetic tendency at the surface of CE charge-ordered manganites 
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Most previous investigations have shown that the surface of a ferromagnetic material may have 
antiferromagnetic tendencies. However, experimentally the opposite effect has been recently ob- 
served: ferromagnetism appears in some nano-sized manganites with a composition such that the 
antiferromagnetic charge-ordered CE state is observed in the bulk. A possible origin is the develop- 
ment of ferromagnetic correlations at the surface of these small systems. To clarify these puzzling 
experimental observations, we have studied the two-orbital double-exchange model near half-doping 
n=0.5, using open boundary conditions to simulate the surface of either bulk or nano-sized man- 
ganites. Considering the enhancement of surface charge density due to a possible AO termination 
(A = trivalent /divalent ion composite, O = oxygen), an unexpected surface phase-separated state 
emerges when the model is studied using Monte Carlo techniques on small clusters. This tendency 
suppresses the CE charge ordering and produces a weak ferromagnetic signal that could explain the 
experimental observations. 

PACS numbers: 75.70.Rf, 75.47.Lx, 75.75.+a 



I. INTRODUCTION 

Perovskite manganites, with a general formula 
AM11O3 , where A is the composite of trivalent rare-earth 
elements and divalent alkaline-earth elements, have at- 
tracted considerable attention since the discovery of the 
colossal magnetoresistance (CMR) effect i Both experi- 
mental and theoretical studies in the last decade have 
unveiled a plethora of phases in manganite compounds, 
with very different macroscopic properties but very sim- 
ilar energies i^i^ The CMR and colossal electroresis- 
tance (CER) 6 effects, which correspond to the obser- 
vation of drastic nonlinear responses of the manganites 
to external stimulations, can be understood as a re- 
sult of the intense competition between the ferromag- 
netic (FM) metallic phase and antiferromagnetic (AFM) 
charge-ordered (CO) insulating phases i ' 8 ' 9 ' 10 

This phase competition is not only sensitive to applied 
external fields but also to the geometric and chemical 
environments of the surface of the system under study. 
In nano-sized materials there is a high surface to volume 
ratio and, as a consequence, the surface effects play a 
crucial role. This influence of the surface has been ob- 
served in a series of recent experiments. On one hand, in 
materials where a FM state is stabilized in the bulk, an 
AFM or spin-glass surface state is foundj 11 ' 12 ! 13 ! 14 ! 15 On 
the other hand, FM tendencies at the surface of nano- 
sized manganites presenting AFM/CO bulk order have 
also been observedffi^^^^^ 3 .^^* 2 ^ 2 ! The first 
tendency toward surface AFM ordering usually appears 
in nano-sized FM or ferrimagnetics materials^ and can 
be understood within the following naive picture. Gener- 
ally, in strongly correlated electron materials, the charge 
conducting properties are determined by the ratio U /W, 
where U is the Hubbard Coulomb repulsion and W is 
the electron hopping bandwidth. At the surface, the 



bandwidth W will be suppressed due to the reduction 
in dimensionality, while the on-site U will not change. 
Therefore, the enhancement of U /W at the surface could 
prefer an insulating surface over a metallic one. The 
insulating phase in manganites is usually AFM. How- 
ever, finding a model and rational for the other ten- 
dency found experimentally, namely a FM tendency at 
the surfaces of CE manganites, is not straightforward. 
A theoretical model that can explain the FM tendency 
at the surface should be able to consider the following 
four experimental signatures. First, the CO phase is 
significantly weakened. The CO transition peak in the 
magnetization (M) vs. temperature (T) curve is sup- 
pressed until it completely disappears with decreasing 
the size of the manganite systems i 16 i 17 i 18 i 19 i 20 i 21 Second, 
M is enhanced at low T and a FM-type hysteresis loop 
is observedi 18 ! 19 !^! 2 !^ 2 .^^! Third, the exchange bias ef- 
fect emerges, indicating a coupling between phases with 
different spin orders i 18 ' 23 And finally, the measurement 
of magnetocaloric properties and low-T specific heat also 
suggest the existence of a FM contribution i 25 ' 26 i 27 

Theoretically, previous investigations have mainly paid 
attention to the surface effects of FM manganite a 29 ' 30 ! 31 
whereas the study of the surface of AFM manganites was 
rare. To understand these phenomena, it is reasonable 
to partition the nano-sized (with typical scales 10 ~ 10 2 
nm) system into two regions: an inner core and a surface 
shell. The physical properties of the inner core should be 
comparable to those of the bulk material. In contrast, 
the properties of the surface shell can be different from 
those of the bulk. The lower coordination number at the 
surface, which effectively reduces the superexchange cou- 
pling, may give rise to a FM tendency at that surface. 
Following this idea, in previous investigations a core-shell 
model was proposed assuming an AFM core wrapped by 
a fully FM surface shell. 32 This model fits some exper- 
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imental results well, but in this scenario it is already 
established from the start the nature of the phases at 
both the core and surface. As a consequence, it is too 
phenomcnological to better understand the true physi- 
cal origin of the AFM-FM transition at the surface. In 
particular, the assumption of a fully FM shell conflicts 
with some experimental results: due to the weak value 
of the magnetization M, the calculated thickness of the 
FM shell can be even thinner than one "molecular" layer, 
which indicates that assuming a fully developed FM spin 
order at the surface is not correct. Thus, theoretical stud- 
ies using realistic microscopic Hamiltonians and unbiased 
assumptions about the surface are necessary to clarify the 
experimental observations found at the surface of nano- 
sized AFM/CO manganites. 

In this paper, the core-shell model is incorporated into 
a two-orbital Hamiltonian for manganites. We show that 
the unscreened Coulomb interactions lead to an increase 
of electron density on the surface. Monte Carlo simula- 
tions reveal that this increase in the density drives the 
surface layer from an AFM/CO state to an unexpected 
phase-separated state, as opposed to a fully developed 
FM state. This surface phase-separated state exhibits 
clear FM signatures, but they are weak, compatible with 
the experimental observations. 



II. MODEL 

To better understand the physics at the surface, here 
we consider a two-orbital model Hamiltonian for mangan- 
ites that includes both finite superexchange coupling and 
the effect of Jahn- Teller phonons. As a well-accepted ap- 
proximation for manganite models, we consider the limit 
of infinite Hund coupling. The Hamiltonian reads: 

a/3 

H = - J2 tfn ijC lc j0 + J AF ]T Si ■ Si 

<ij> <ij> 

+ p) n i + A /XQim + Q2iT xl + Q3iT zi ) 

i i 

+^E( 2 ^ + ^ + Q^)- (i) 

i 

Here, the first term is the two-orbital double exchange 
interaction, a and (3 denote the two Mn e g -orbitals a 
(d x 2_ y 2) and b (d3 2 2_ r 2). c,; a (cj a ) annihilates (creates) 
an e g electron in orbital a of site i with its spin paral- 
lel to the localized t 2g spin Si. The hopping direction is 
denoted by r. As discussed in previous literature^ the 
hopping amplitudes are t a x a = t a y a = 3t b J> = 3t bb = 3/4, 

tf = t b y a = -tf = -4° = v/3/4, t a z a = tf = t h z a = 
and t b z b — 1 (energy unit). The infinite Hund cou- 
pling generates the factor fiy = cos(6>j/2) cos(#j/2) + 
sin(6>j/2) sin(#j/2) exp[— i(^ — <j>j)], where 8 and are 
the angles of the t 2g spins in spherical coordinates. The 
second term is the superexchange interaction between 
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FIG. 1: (Color online) Sketch of the model and geometry 
used in this investigation. Left: a cubic lattice with one open 
surface (yellow layer). All spins in the bottom layer (orange) 
are frozen to the CE-type AFM pattern. Right: the chemical 
unit taken from the surface layer of the cube at the left. The 
AO sheet termination is considered here. Thus, the surface 
formula is Ai.^MnOa.^. 

nearest-neighbor (NN) t 2g spins. In the third term, fi 
is the chemical potential, q corresponds to a site depen- 
dent Coulomb potential. The origin and relevance of this 
term will be discussed in detail later in this section, rii is 
the e g charge-density at site i. The fourth term stands for 
the electron-phonon coupling. The Qs are phonons cor- 
responding to Jahn- Teller modes (Q 2 and Q3) and the 
breathing mode (Qi). r is the orbital pseudospin opera- 
tor, giving t x = c^Cb+c^Ca and t 2 = c> a Ca-c\c b . The last 
term is the elastic energy of the phonons. For simplicity, 
we have already assumed in this model that both the tig 
spins and the phononic degrees of freedom are classical 
variables. The above described Hamiltonian is solved via 
a combination of exact diagonalization and Monte Carlo 
(MC) techniques: classical ti g spins and phonons evolve 
following the MC procedure; and at each MC step, the 
fermionic sector of the Hamiltonian is numerically ex- 
actly diagonalized. The first 10 4 MC steps are used for 
thermal equilibrium and another 6 x 10 3 MC steps are 
used for measurement. More details of the Hamiltonian 
and MC technique can be found in Refs. 2,3. 

Since a FM tendency is often found in half- 
dope d 16 i 17 ' 18 i 24 or nearly half-dope d 19 ' 26 nano-sized man- 
ganites, and since a CE-type AFM/CO phase usually ap- 
pears in half-doped narrow band manganites)^ 3 . we solve 
the two-orbital model using densities corresponding to 
half-doped systems and for couplings where a bulk CE 
phase exists. This can be done by tuning the couplings 
(Jaf, A) to be within the CE regime of the phase diagram 
obtained in previous studies . 34 ' 35 Although only one set 
(Jaf, A) will be investigated in the studies described be- 
low, we believe our qualitative conclusions remain valid 
for other choices of parameter sets, as long as they are 
within the the CE region of the phase diagram. 

Although there are several possible surface directions 
in real cases, here only the simplest case, i.e., the (001) 
surface, will be considered, similarly as in most for- 
mer theoretical investigations in this context. 29 ' 30 ' 31 ' 36 
Strictly speaking, the surface problem should be consid- 
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ered on a three-dimensional half-infinite cubic lattice, e.g. 
infinite in the x and y directions, but semi-infinite in the 
z direction. In practice, to address this problem numeri- 
cally the above two-orbital model Hamiltonian is studied 
using a L x L x L cubic lattice, as shown in Fig. 1. Peri- 
odic boundary conditions (PBC) are applied in both the 
x and y directions. Following the core-shell-model idea 
discussed in the previous section, in the z direction the 
spins in the bottom layer (z = L) are here fixed so that 
they have the same CE-type AFM pattern as in the bulk. 
By this procedure, the properties of the half-infinite core 
are encoded in the bottom layer in our model. Hence, 
the rest L — 1 layers are to be considered as the outer 
shell. To better justify the above assumption, L must be 
large enough that surface effects are limited within the 
outer L — 1 shell layers. Our results below suggest that 
in practice L — 4 is enough for these purposes. This is 
fortunate, since the numerical studies described here are 
rather CPU time consuming. 

For the surface layer (z = 1), we have to take into 
account the "termination" procedure. In this paper, a 
clean-limit AO sheet is considered as the termination, 
which makes the "molecular composition" of the sur- 
face layer to be A1.5MnO3.5- Therefore, the outmost 
AO sheet transfers an extra 0.25 electron per site to the 
nearby Mn cations, and it is positively polarized. The 
polarized AO sheet introduces an unscreened Coulomb 
attraction on the surface. As a result, those extra elec- 
trons must be localized near the surface. The effect of 
this surface Coulomb attraction has been taken into ac- 
count here via an effective negative potential near the 
surface: ej = V < when i belongs to the surface layer 
and ej = otherwise. This may be justified because 
of the typical short Thomas-Fermi-like screening length 
found in many of these materials. The boundary con- 
dition on the surface layer is then set to be open. The 
consideration of the AO sheet has another important ef- 
fect: it keeps the oxygen octahedrons complete for the 
outmost Mn cations. Therefore, the phonon modes do 
not need to be changed, even for the surface layer. For 
the phonons, PBC are used in all directions for simplicity. 
This will not affect much the central physics discussed in 
this manuscript since the oxygen displacements along the 
z axis are negligible in both the FM and CE type AFM 
phases^ In addition to the AO termination, the other 
possible choice for the surface is a clean-limit Mn02 sheet 
as the termination, but this will not be considered in our 
current model. However, the possible effects of this al- 
ternative termination will also be discussed in the next 
section. 

To contrast results with those of the open surface case, 
the model will also be studied following assumptions that 
address the bulk material: this corresponds to using PBC 
in all directions, for both spins and phonons, and with- 
out freezing spins anywhere. Both simulations are per- 
formed at T = 0.02 (experimentally, it corresponds to 
50 ~ 100 K). To simulate the bulk material, we set 
e, = for all sites. The parameters Jaf — 0.1 and 



A = 1.2 are used to obtain a stable CE phase with av- 
erage density < rii >= 0.5. Note that it is well known 
that the CE phase is stable over a broad range of cou- 
plings of the half-doped manganites, thus this selection 
of parameters should not be considered arbitrary or fine 
tuned. The stability of the CE phase is also confirmed in 
the present simulations by analyzing the spin and charge 
structure factors, which will be discussed in detail in 
the next section. As observed in Fig. 2(a), this corre- 
sponds to choosing the chemical potential in the win- 
dow — 1.1 < /i < —0.95. These parameters will also be 
adopted in our subsequent simulation of the open surface 
model. 

To obtain reasonable results from the simulation of the 
open surface model, we have to set the Coulomb potential 
to an appropriate value. Actually, an optimal value of V 
exists in order to fulfill the following three criteria: first, 
the average density per site of the entire Lx Lx L system 
< n, > should be equal to (0.5 + 0.25/L) = 0.5625 to 
keep the charge neutral; second, the chemical potential 
should remain within (—1.1, —0.95), which is required to 
have the CE phase stable far from the surface; third, the 
charge density in the bottom layer should be very close 
to 0.50 to match the frozen CE type. In order to find 
the optimal V, we tested several values from to —0.6, 
stepped by —0.1. In this range, V = —0.4 was found to 
be a proper parameter at /1 = —1. Thus, V = —0.4 was 
the value adopted in the following simulations. 



III. RESULTS AND DISCUSSION 

Using as initial configuration a perfect CE spin 
pattern^ and random phonons, the MC simulation was 
carried out under the assumption V = —0.4 for the open- 
surface model. The averaged e g charge-density as a func- 
tion of /i is shown in Fig. 2(b), allowing us to compare the 
results for the bulk material model with those of the open 
surface. An interesting result is that the n=0.5 plateau 
in the bulk model, which corresponds to the stable CE 
phase, disappears in the open-surface model. The charge 
density increases with increasing chemical potential but 
it has large fluctuations (see error bars), suggesting that 
the system presents a strong competition between phases 
with different densities. In the following, we will focus on 
the properties at \x = —1, with the average total charge 
density close to 0.5625^ 

To understand the origin of the different average charge 
density between bulk and open-surface models shown in 
Fig. 2(b), it is important to analyze the e g charge-density 
at each layer. The results are shown in Fig. 3(a). The 
charge density is almost exactly 0.50 at the bottom layer 
and fluctuates around 0.50 in the two middle layers. The 
most prominent change occurs at the surface layer, where 
the charge density increases to 0.75, due to the presence 
of the Coulombic term. As expected, the extra 0.25 elec- 
tron from the outmost AO sheet is mainly located in 
the first Mn sheet, which offsets the Coulomb interaction 
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arising from the outmost ^40 sheet for the Mn cations 
of the second layer. This is consistent with ab-initio 
calculations that show that for the ^40 termination the 
uncompensated electrons are accumulated mainly at the 
surface.- 5 This result is also consistent with the assump- 
tion that €i is nonzero only for the first layer. 

It should be noted that the density 0.75 usually corre- 
sponds to a FM phase in bulk manganites, as indicated 
in Fig. 2(a). If the bulk phase diagram remains valid at 
the surface, we would expect a FM state there. Then, 
the FM tendency could be naturally explained as a re- 
sult of a density-driven transition. To verify this possi- 
ble FM tendency, the average NN spin-correlation as a 
function of layer index is shown in Fig. 3(b). Starting 
from the bottom layer (L = 4), we find that the in-plane 
NN spin correlations are almost zero for layer L = 4 and 
L = 3, implying that the FM and AFM links have almost 
the same population. The inter-layer correlation between 
these two layers is close to —1, indicating a fully AFM 
connection between the two layers. These are consistent 
with the picture that a CE AFM state is stabilized in 
these two layers, given that the charge density is about 
0.50. Interestingly, we do see an increase of both the 
in-plane and inter-layer NN spin-correlations as we ap- 
proach the surface. Both the correlations in the first layer 
and between the first and second layers take positive val- 
ues, and they display a clear FM tendency at the surface 
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FIG. 3: (a) The e g charge-density for each layer (counted 
from the surface), (b) The corresponding averaged NN spin- 
correlation, within (integer index) and between (half-integer 
index) layers. 
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FIG. 2: (Color online) The total e g charge-density as a func- 
tion of the chemical potential using L = 4. (a) is for the 
case of the bulk material model. The corresponding NN 
spin-correlation is also shown. There is a < rn >= 0.50 
plateau with the NN spin-correlation ~ —0.3 in the range 
— 1.1 < fl < —0.95, indicating a stable CE phase, (b) con- 
trasts the results for the open-surface and the bulk material 
models near half-doping. 



layer. But the positive value for the in-plane correlation 
is rather small (~ 0.15) at the surface layer, which sug- 
gests that the state is only partially FM. Therefore, the 
idea of a density-driven transition is too simplistic and 
not quite correct. This already shows an interesting con- 
clusion of our research: the phase diagram at the surface 
cannot be obtained by merely analyzing the bulk phase 
diagram at the appropriate charge density but a special 
investigation is needed to clarify the surface's properties. 

To visualize the nature of the surface state more ex- 
plicitly, we study the distribution of local charge density 
on the surface layer. In a FM phase, the distribution of 
local charge density n, is approximately uniform, in con- 
trast to the charge disproportion typical of an AFM/CO 
phased The distribution of local charge density on the 
surface layer of the model studied here is presented in 
Fig. 4(a), and the result at the bottom layer is also shown 
in the same figure for contrast. From the regular charge 
pattern, it is clear that the CE charge ordering in the 
bottom layer is very stable. Since the CE spin pattern 
is fixed at the bottom layer, the charge disproportion 
between large-density sites (bridge-sites of the zig-zag 
chains) and smalle-density sites (corner-sites of the zig- 
zag chains) will not be smeared by MC average. However, 
for the surface layer the charge distribution is not uni- 
form, and it is not regularly distributed. In some sites, 
the densities are large and close to 1, but in other sites 
the densities can be as low as approximately 0.5. This 
inhomogeneous distribution persists prominently even in 
the MC-averaged result, which rules out the possibility 
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of observing a non-uniform charge state due to thermal 
fluctuations. This inhomogeneity can also be confirmed 
from the orbital occupation, as shown in Fig. 4(b). In 
contrast to the distinct CE type orbital-ordering in the 
bottom layer, the orbital distribution in the surface layer 
shows two regions corresponding to charge inhomogene- 
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FIG. 4: (Color online) (a) The e 3 charge-distribution corre- 
sponding to the surface layer (left) and bottom layer (right). 
Here are shown both a typical MC "snapshot" (upper panels) 
and MC averaged results (lower panel, averaged over 6, 000 
MC steps). The size of the circles is in proportion to the 
local charge density m. For the surface layer, the dotted 
lines separate the high density (n, > 0.75) and low density 
(rii < 0.75) regions, (b) Sketch of the orbital occupation 
(based on the MC-averaged r X i and r z i) corresponding to the 
above-described MC-averaged charge distribution, (c) MC- 
time evolution after thermal equilibrium of the e 9 charge- 
density (average value for all the sites). The frequent tunnel- 
ing events are prominent among the possible densities, indi- 
cating tendencies toward electronic phase separation, which 
is a first-order transition when fi is varied. 



ity. The high-density region shows an orbital ordering 
similar to the case in undoped manganites, in contrast to 
the low-density region which shows the orbital disorder 
similar to the case in FM manganites. In addition, the 
MC-time evolution of the e g charge-density is presented 
in Fig. 4(c), showing that the tunneling events are promi- 
nent among several possible densities. These tunneling 
events are characteristic of a first-order phase transition 
varying /i, rather than standard thermal fluctuations. 41 
Therefore, the inhomogeneous charge distribution at the 
surface should be attributed to tendencies in the model 
toward nanoscale electronic phase separation;^ similarly 
as those observed in bulk simulations in other regions of 
parameter space. 

To quantitatively reveal the competing phases that ap- 
pear at the surface, the spin structure factor S'(q) for 
each layer is calculated via a Fourier transformation of 
the spin correlation function) 4 ^ 



L 4 



-iv) 



(2) 



By monitoring the q dependent spin structure factor, we 
may detect possible weak spin-order signals in the com- 
plicated real-space spin patter n 35 i 43 because each spin 
order corresponds to a unique set of characteristic q 
vectors.— In Fig. 5(a), the S(q)s for some possible com- 
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FIG. 5: (Color online) (a) The spin structure factor values 
for several spin orders, showing the reduction in the tendency 
toward a CE state as the surface is reached. At this sur- 
face, the FM and E phase tendencies are dominant, (b) The 
charge structure factor at (7r,7r) as a function of the layer in- 
dex, showing the reduction of the CE-phase staggered charge- 
order tendencies as the surface is reached. 
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ponents (summed over all characteristic q vectors for 
each spin order) are shown layer by layer. For the bot- 
tom layer, which is frozen into a CE spin pattern, both 
the C and E components contribute 50%, respectively, 
as expected. The CE phase tendencies become gradually 
weaker with decreasing layer index, indicated by the de- 
creasing values and large fluctuations of the correspond- 
ing 5(q)s. At the surface, the C component is very weak 
(~ 10%) and the dominating components here are FM 
(~ 40%) and E (~ 30%). In this case, the E value does 
not match the C value anymore because the site regions 
where rij ~ 1 can also contribute to the E type order 4^ 
In addition to the S(q) spin structure factor, the charge 
structure factor C(q) in each layer is also calculated to 
characterize the charge ordering: 

= Jji X> - ~ n0e <q - (p,-ri) , (3) 

where ni is the average density of each layer. Here, 
only C(iy,tt) is shown in Fig. 5(b) since all other com- 
ponents are very close to zero. For the CE phase at the 
bottom layer, C(ir,ir) is about 18%, a result consistent 
with the expected charge-ordering with charge dispro- 
portion ~ 0.4 (the charge difference between high- and 
low-density sites). Decreasing the layer index, C(7r,7r) 
decreases monotonically as we move toward the surface, 
until it completely disappears at the surface. This is 
also straightforward to understand since charge-ordering 
usually accompanies the AFM phases instead of the FM 
ones, even in bulk manganites. 

In the above simulation, the ground state is a robust 
CO CE phase by fixing the proper Jaf and A. There- 
fore, it is worth to address the other possible cases in real 
half-doped manganites. Here we will give a brief analysis 
based on the half-doped phase diagram obtained in previ- 
ous works i 34 ' 35 On one hand, by decreasing Jaf only, the 
ground phase can change from the CE to A-type AFM, 
then finally to FM phase. This process corresponds to 
the experimental observed transition from the narrow 
band manganites to middle band one, then finally to 
wide band manganties4^ On the other hand, by decreas- 
ing A only, the charge-ordering (or the degree of charge 
disproportionation) will become weaker and weaker, till 
completely turn to the FM phase when (Jaf, A) crosses 
the phase boundary between CE and FM phases. In 
both cases, FM tendency will be enhanced. In short, 
our above simulation mainly aims at the family of nar- 
row band manganties whose CO CE phase is stable, e.g. 
Nd .5Cao.5Mn0 3 and Pr .5Cao.5Mn0 3 . 

In the simulations described above, only the AO sheet 
termination was considered. But it is necessary to discuss 
the other choice already mentioned, namely a Mn02- 
sheet termination. In the case of this Mn02 termina- 
tion, the cubic symmetry is lost at the surface due to 
the breakdown of the oxygen octahedrons. Therefore, 
the 3d energy levels of the Mn cations at the surface are 
different from those with full oxygen octahedrons. In 



particular, the energy of the d3 Z 2_ r 2 orbital will be low- 
ered substantially. An early model study by Calderon 
et al showed that the Mn02 termination would gener- 
ate an AFM surface for FM manganites because the e g 
density at the surface was enhanced from ~ 0.7 to l3i 
However, an ab-initio calculation by Fang et al. showed 
the reverse result: decreased e g density at the surface 
by Mn02 termination^ Therefore, whether the Mn02 
termination can generate the FM tendency, e.g. by en- 
hancing the surface charge density from 0.5 to 0.75 as it 
occurs in the case of the AO termination, remains unclear 
and is an interesting subject of investigations. However, 
this M11O2 termination is far more complex than the case 
studied here, and beyond the scope of the present work. 

It is important to remark that some other extrinsic 
factors, such as defects of cations/oxygen and recompo- 
sition of surface structures, may also affect the physical 
properties of real manganites. Even qualitatively con- 
sidering these effects, our model still gives reasonable 
results. For a crude comparison with experiments, the 
FM component fraction can be estimated as 40% of the 
surface layer (Fig. 5(a)). Therefore, to compare with ex- 
periments the FM fraction predicted by our study we 
should use the number 40% x surface/ volume (at zero 
magnetic field and low T), where surface/ volume should 
be calculated based on the actual shape and size of the 
nano-sized clusters used experimentally. By this pro- 
cedure, our theory agrees with the weak magnetization 
found in Ndo.sCao.sMnOs nanoparticles (diameter ~ 20 
nm, experimental magnetization M ~ 3 emu/g ~ 3% 
of the saturation magnetization M s , while in our model 
our estimation gives ~ 5%M s )i£ and Pro.sCao.sMnOa 
nanowires (diameter ~ 50 nm, experimental magnetiza- 
tion M ~ 1 emu/g ~ 1%M S , while our model estimation 
is ~ 1%M S )M It should be noted that the phenomeno- 
logical core-shell model can not explain these very weak 
magnetizations^ The agreement between our estimates 
and experimental data suggests that our model at least 
grasps the main physics of the surface effects in the nano- 
sized CE-phase manganites. 



IV. CONCLUSION 

In conclusion, we have performed a Monte Carlo study 
of the CE-type AFM/CO phase in a 3D lattice with 
an open surface. The AO sheet termination leads to 
the generation of an extra 0.25 electron-per-site at the 
surface, here simulated by the introduction of an un- 
screened Coulomb attraction at that surface. As a re- 
sult, the charge density on the surface Mn-layer was en- 
hanced from 0.50 to ~ 0.75. The charge density ~ 0.75 
usually corresponds to a fully FM phase in bulk man- 
ganites. However, within the Monte Carlo simulations 
for small clusters discussed in this manuscript, the sur- 
face was found to have a nontrivial nanoscale electronic 
phase separated state. At the surface, the charge dis- 
tribution was found to be inhomogeneous and coexist- 
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ing with a weak FM spin correlation. The studies of 
both the charge structure factor and spin structure fac- 
tor confirmed the suppression of AFM/CO order and the 
enhancement of FM order near the surface. Our result is 
helpful to understand the weak FM tendencies observed 
in nano-sized AFM/CO manganites. However, clearly 
these results have to be considered as just a first step 
toward the understanding of the phase diagram at the 
surfaces of manganites. Larger clusters and other nu- 
merical/analytical techniques should be used to confirm 
our results and further explore the physics unveiled here. 
Also, a systematic study varying parameters of the many 
tendencies expected in the anticipated rich phase dia- 
gram of these compounds at the surface should be carried 
out in future investigations. 
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